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Abstract 

The contribution of the tensor interaction in Hartree-Fock calculation for closed 
shell nuclei is studied. The investigation is done neglecting the spin-orbit component 
of the force. As expected from nuclear matter estimates we find that the effect of 
the tensor interaction is negligible in the case of spin-isospin saturated nuclei and 
also for medium and heavy nuclei with closed shells. In light nuclei with only one of 
the two spin-orbit partners filled, such as 12 C and C, the tensor interaction plays 
an important role. 
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1 Introduction 



One of the peculiarities of the nucleon-nucleon interaction is the presence of non central 
terms. Already at the beginning of the 40 's the existence of a quadrupole moment of the 
deuteron was recognized This fact was explained by Rarita and Schwinger in terms 
of a static tensor force ||. The presence of tensor terms is an essential ingredient in the 
modern phase-shift analysis of nucleon-nucleon scattering data || . In a meson exchange 
picture of the nucleon-nucleon interaction, tensor terms arise already when the exchange 
of the lightest and better known meson, the pion, is considered. 

The relevance of the tensor part of the nucleon-nucleon interaction is not only restricted 
to two-body systems. Calculations based on microscopic nucleon-nucleon interactions of 
light and medium-heavy nuclei Q and nuclear matter || show that binding is obtained 
only because of the presence of the tensor force. 

Although the tensor terms of the interactions are essential in microscopic calculations, 
they are usually neglected in effective theories where the Schrodinger equation is solved 
at the mean-field level, since their contribution is though to be small. This attitude is not 
unreasonable if one consider that the tensor terms of the interactions do not contribute in 
mean-field calculations of the binding energy of symmetric nuclear matter. On the other 
hand, since the spin and isospin structure of finite nuclei is richer than that of nuclear 
matter, one may expect to find, in the former systems some observable sensitive to the 
presence of the tensor terms of the interaction. In this respect Random Phase Approx- 
imation calculations have shown that magnetic excitations can be reasonably described 
only when tensor terms are included in the effective interaction ||. 

We have investigated the effects of the tensor terms of the nucleon-nucleon interaction 
on ground state properties of doubly-magic nuclei and in this work we would like to answer 
the question whether it is reasonable to neglect this part of the interaction in Hartree- 
Fock calculations. We started our study by using a well know finite range interaction, 
the Brink and Boeker Bl force J7|, defined only in the four central components. To 
this basic interaction we have added the tensor and tensor-isospin terms taken from the 
realistic Argonne Vu potential || and from the phenomenological Jiilich-Stony Brook 
(JSB) interaction ||. Our investigation consists in comparing the results obtained with 
these three different interactions for various closed shell nuclei. 

We briefly describe the formalism in the next section and we present the results of 
our investigation in Sect. |3[ We shall show and explain the fact that the contribution of 
the tensor interaction is negligible for those nuclei where all the spin-orbit partners single 
particle levels are occupied. For the other nuclei, these contributions remains small if the 
nuclei are heavy. 

2 The formalism 

The Hartree-Fock (HF) equations are obtained searching for a minimum of the energy 
functional within a Hilbert space restricted to states which are Slater determinant of single 
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particle (sp) wave functions. In this way the many-body problem is transformed in a set 
of many one-body problems. For hamiltonians containing only two-body interactions, 
without explicit density dependent terms, one has to solve, in coordinate space, a set of 
equations of the type: 

-7^Vi 2 fc ( ri ) + C/(n)0 fc (n) - / d 3 r 2 ^(n,r 2 )(/> fc (r 2 ) = e fc fe (n), (1) 
zrrik J 

where we have indicated with k the set of quantum numbers characterizing the sp wave 
function (fik and the sp energy e&. The two quantities U and W are the standard Hartree 
and Fock-Dirac potentials || 

The only input of the theory is the effective nucleon-nucleon interaction which we 
have chosen cLS cL SIX components finite-range two-body force of the type: 

V(r u r 2 ) = £ y p ( ri ,r 2 )O p (l,2) (2) 
P =i 

where O p (l,2) indicates the following operators 1, T\ ■ r 2 , cr x • <x 2 , <j x ■ <j 2 T\ ■ r 2 , S\ 2y 
S12T1 ■ r 2 , for p = 1, ... ,6, respectively. In this expression we have adopted the usual 
convention for the tensor operator: 

c a (o"i-ri2)(o-2-ri 2 ) 

&12 = 3 <Tl • <T 2 , (3) 

with r 12 = i"i — r 2 . 

We suppose that all the functions V p in eq. (0) depend only from the relative distance 
between the two particles, |i"i 2 |. It is therefore convenient to use expressions of the sp 
wave functions where the angular and radial coordinates are separated: 

vLi At) 



J2(l^s\jm)Y ¥ (e,(f)xsXt- (4) 

/is ^ 



In the above expression 9 and (p are the angular coordinates, is a spherical harmonics, 
Xs and Xt are the spin and isospin wave functions, s and t represent the third components 
of spin and isospin, and (lp^s\jm) is a Clebsh-Gordan coefficient. 

Multiplying eq. ([I]) by (Ikjk^ktkl and integrating on the angular, spin and isospin 
variables we obtain: 



% 2 ( d 2 i k (i k + iy 



U n k l k j k ( r ) + ^n k l k j k ( r ) U n k l k i k ( r ) ^n k l k j k ( r ) 



Irak \dr 2 r 2 ) fe/fc:,fe nkhik" ' n k i k j k v ) ^n k i k j 

tu th- 
c ** If 

-rikhik n k l k j k 



e n k l k i k U n k l k j k ( r ) ■ (5) 



To evaluate U^ k i k j k (r) and (r) we write the interaction 1/(r!,r 2 ) in terms of 

its Fourier transform. This allows us to separate the functional dependence of the two 
coordinates ri and r 2 . 
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Spin and isospin conservation implies that in the direct term only the scalar (p=l) 
and the isospin (p=2) channels of the interaction contribute. For these direct terms we 
find the expression: 



(6) 



KXi^r)}^) = 4v ^ J drV 2 V p (r,r') [U k {r')] {p) 
where we have defined 

V£(r,r') = J dqq 2 j L (qr)j L (qr')Vi°\q), 
with ji{x) a spherical Bessel function and 



(7) 



V}°\q) 

Finally, the function U is: 



dr r 2 jo(qr) V p (r) 



p(r) , p = 1 

2t k [ff(r)-p n (r)} , p = 2 



with t k = | ( — |) according to k being a proton (neutron). 
In eq. (p|) we have indicated with p the nucleon density: 



p( r ) = 



nljt 



2j + l ( <ij(r) \ 2 
4-7T I r I 



(8) 



(9) 



(10) 



and with fp and p n the proton and neutron densities, defined as in eq. ([H]) without the 
sum on t, such as p = fP + p n . 



All the channels of the force contribute to the exchange term in eq. 
expressed as: 



which are 



E 



rikhjk K 



(p) 



£(i k + l i + L) {2 3l + 1) (2L + 1) I p [£ kiL {r, r')\ p) , 



11' 



where we have introduced the function 



(12) 



The term I p is obtained from the calculation of the isospin matrix elements and is given 
by: 



6, 



p = 1,3,5 



2S tk ,. u + 6 tk , u p = 2,4,6 ' 



(13) 
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where we have indicated with 8 the Kronecker symbol. Finally, the function £ includes 
the angular momentum coupling pieces. For the scalar and isospin terms we have the 
following expression: 

2 



3k 3i L 
1 1 



V p L (r,r') 



[Skidr, 0^=1,2) 

V 2 2 

where we have made use of the Racah 3j symbol. 

The contribution of the spin dependent terms is given by: 

[Skidr, r')] (p=3i4 ) 



(14) 



Ik h L 




3k 3i L 

2 2 u 



V£(r,r') 



(15) 



Finally, for the tensor channels we have: 

[^,L(r,r')] (p=5j6) = 2V30(2Z, + 1) (2l k + 1) £ (-i) 1 ' 1 ' 



-l) J V£ L ,(r,r') 



L'J 



(2J+1) (21/ + 1) 





V 






L 




{i 


1 


5) 


{« 





:) 



(16) 




l k li L 




where the Racah 6j and 9j symbols have been used and where we have introduced the 
function 

' Mq 2 3L(qr)ju(qr')V^ 2 \q), (17) 



with 



dr r 2 j 2 (qr) V p {r) 



The solution of eq. @ provide us with the sp energies e and wave functions u(r). 
This solution has been obtained iteratively using the plane waves expansion method of 
refs. [10]. The center of mass motion has been considered in its simplest approximation, 
consisting in inserting the nucleon reduced mass in the hamiltonian. The sp wave functions 
used to start the iterative procedure have been generated by a Saxon-Woods potential 
without spin-orbit and Coulomb terms. This means we have used the same starting wave 
functions for both protons and neutrons and for sp states which are spin-orbit partners. 

The total energy of the system is calculated as: 

1 ^ r dVadVa^CrO^^^drjaD^Cn)^^) 



E 



e k 
k 



E 

ki 



d 3 n d 3 r 2 0*( ri ) 0*(r 2 ) V{\v l2 \) ^(n) M^) 



(19) 



E (2J + 1) e 

nljt 



t 

nl j 



An drr p(r 




E(2j + 1) / dr EK< 

nljt \p=l 



(20) 
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3 Results and discussion 



Our calculations are based on a finite range effective force widely used in HF calculations, 
the Bl interaction of Brink and Boeker which was constructed to reproduce the exper- 
imental binding energies of 4 He, 16 and 40 Ca. With this interaction we have performed 
HF calculations also for 12 C, 14 C, 48 Ca and 208 Pb nuclei. 

A comparison between the binding energies obtained in our calculations with the 
experimental values [1 1] is given in table |l|. While the energies of 4 He, 16 and 40 Ca are 



rather well reproduced, the model fails in describing those of the other nuclei. A study 
of the contribution of the four components of the force shows that only the spin-isospin 
(p=4) term provides attraction. The contribution of the isospin (p=2) and spin (p=3) 
channels is the same for the Z = N nuclei. Since in the present calculations protons 
and neutrons sp wave functions are identical, for Z = N nuclei the proton and neutron 
densities are equal, therefore the direct isospin term, eq. (^), is zero and the two exchange 
terms for p=2 and p=3, give the same result. Of course this picture breaks down for nuclei 
with neutron excess like 14 C, 48 Ca and 208 Pb. 

To investigate the influence of the tensor component of the force on HF calculations 
we have added to the Bl interaction the tensor force obtained from the Argonne 
potential || and that obtained from the Jiilich-Stony Brook (JSB) interaction ||. The 
tensor terms of the JSB interaction are constructed by considering the exchange of the ir 
plus p mesons, therefore they are active only in the tensor isospin channel. On the other 
hand, the pure tensor term of the Argonne potential is very small if compared with the 
tensor-isospin one. For this reason in our calculations we have used only the tensor-isospin 
terms of the two interactions. These terms are shown in the panel (a) of fig. [I]. 

This treatment of the tensor part of the interaction reproduces rather well the tensor 
terms of an effective interaction calculated within a Brueckner G-matrix approach. In ref. 
[ |T2| is it shown that the effective force in the tensor channels is very similar to the bare 
nucleon-nucleon interaction. From fig. p] the similarity between the JSB and Argonne 
tensor isospin channels is evident in the low momentum region. 

We did not use the full Argonne and JSB interactions because they are not suited for 
HF calculations. The Argonne potential is a realistic nucleon-nucleon interaction which 
reproduces the experimental data of the two nucleon systems, and therefore it has a strong 
repulsive core. With this interaction we have calculated 16 and 40 Ca and both of them 
are unbound (52.08 and 217.18 MeV respectively). The JSB interaction is an effective 
force designed to describe low-lying excited states within the Landau-Migdal theory of 
finite Fermi systems, and for HF calculation it turns out to be too attractive. 

The binding energies obtained with the interactions constructed by adding to the 
Bl interaction the tensor terms of the Argonne Vu potential (A) and that of the JSB 
interaction (JSB) are shown in table [|. A comparison with the results of table [l] shows 
that the contribution of the tensor force is very small. 

The set of sp energies for the nuclei we have studied is shown in tables |3|-[?J For brevity 
we do not show these energies for the 208 Pb case. One can notice that in the 16 and 
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40 Ca nuclei, the sp energies of the spin-orbit partner levels are the same, while for 14 C 
and 48 Ca they are different. 



This results can be understood analyzing eq. (0). If the r-dependent terms would be 
constant, the angular momentum sums for nuclei having spin-orbit partner states occu- 
pied would be exactly zero. In reality each term is identified by small variations produced 
by the integral on q, depending from L and V and the integral on r depending from 
the sp wave functions. These variations are not large enough to prevent big cancelations 
between the contribution of the spin-orbit partners terms. The contribution of the tensor 
interaction can be significant in nuclei where not all the sp spin-orbit partners levels are 
occupied. In closed shell nuclei there are at most two of such a levels, one for neutrons and 
one for protons. The contribution of the tensor channel should be compared with that of 
the other channels where all the sp states contribute. In light nuclei, like 12 C, the tensor 
interaction could produce a noticeable contribution if compared with that generated by 
the other channels of the force which involves few sp levels. We expect that in heavier 
nuclei the tensor contribution produced by only one sp level should be compared with 
that produced by the other, many, sp levels for the central channels of the force. For this 
reason in the tables the effect of the tensor force is noticeable only for those nuclei 
having spin-orbit partners not fully occupied. 

It is interesting to notice that the spin-orbit splitting produced by the tensor force 
goes in the opposite direction with respect to the experimental one. Clearly, a more 
realistic description of the nuclear ground state in terms of mean field model requires the 
presence of explicit spin-orbit terms in the interaction. 

In order to test the sensitivity of the 12 C ground state properties to the tensor inter- 
action we have inserted a tensor-isospin channel of gaussian shape (in coordinate space) 
whose parameters have been fixed to reproduce the experimental binding energy of this 
nucleus. This new tensor force is shown in the panel (b) of fig. |I| and the results obtained 
with this new interaction for the other nuclei are shown in the tables by the columns 
labeled T6. It is worth to notice that this interaction has opposite sign with respect to 
the tensor terms of the JSB and Argonne interaction and a much stronger strength. 

The comparison done in table ^| between the binding energies obtained with this new 
interaction with those previously obtained shows that only the results of the two carbon 
isotopes have been strongly modified. In this table we see that the binding energies of the 
nuclei with all spin-orbit partners levels occupied are very little modified by the tensor 
force. The 48 Ca and 208 Pb nuclei show more sensitivity to the modification of the tensor 
component of the force, but it is however rather small. The same kind of effects can be 
seen in the other tables showing the sp energies. The sp energies of 12 C and 14 C nuclei are 
strongly modified by the new tensor force, but those of 16 and 40 Ca remain the same. 
In 48 Ca the modification of the sp energies are relevant, even though the change on the 
binding energy is very small (~ 2%). 

In order to see how the modification of the tensor interaction can affect the sp wave 
functions we have calculated the nucleon distributions. As example of the results obtained 
we show in fig. the proton distributions of the 12 C, 14 C, 40 Ca and 48 Ca nuclei. On the 
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scale of the figure the curves representing the Argonne and JSB results are overlapped 
to the curves of the Bl results. The figure shows that the JSB and Argonne tensor 
potentials do not produce sizeable effects on the proton distributions. The situation is 
quite different when the T6 tensor terms is used. The proton distribution of 12 C and 14 C 
are quite different, while the effect on 40 Ca and 48 Ca are negligible. These results confirm 
the previous analysis on the behaviour of the tensor interaction. The tensor force plays 
an important role in 12 C and 14 C, but its effects are essentially zero for nuclei with fully 
occupied spin-orbit partner levels, and very small for medium-heavy nuclei. 



4 Conclusions 

In this work we have studied the effects of the tensor interaction on HF calculations. We 
found that there are big cancelations between the contributions of sp levels which are 
spin-orbit partners. For this reason, only the contribution of those sp states not having 
occupied the spin orbit partner level is significant. A closed shell nucleus can have at 
most two of these states, one for protons and one for neutrons, and, as a consequence, 
the effect is relatively small in heavy nuclei where it has to be compared with the effect 
of the other terms of the force active on all the sp levels. The situation is different for 12 C 
where only two levels are occupied, therefore, as we have shown, the effect of the tensor 
interaction is comparable with that of the other terms. 

One may argue whether the difficulties in reproducing the 12 C properties of traditional 
HF calculations can be overcome by inserting the tensor terms of the interaction. We have 
shown that, if these terms are included, it is possible to reproduce the 12 C binding energies 
without spoiling the agreement obtained for the other nuclei. To obtain this results the 
use of tensor forces very different from those of Argonne and JSB interactions has been 



necessary. On the other hand, the study of excited low-lying magnetic states gives 
strong indications that the effective tensor isospin interaction should be quite similar to 
those of the Argonne and JSB interactions. The unified description of ground and ex- 
cited nuclear states within a Hartree-Fock plus Random Phase Approximation theoretical 
scheme, rules out the possibility of fixing the tensor part of the effective interaction with 
a fit to the 12 C ground state properties. 

The task of constructing a new finite-range interaction to be used in mean field cal- 
culations was beyond the aim of the present paper, and it will be pursued in the future, 
by adding finite range spin-orbit and density dependent terms. In any case the results 
we have presented show the need of including tensor terms in this effective interaction if 
one wants the model to be valid on the full isotope table. 
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Table 1: Binding energies, in MeV, obtained with the Bl interaction. The experimental 
values are taken from ref. [11]. 





exp 


Bl 


4 He 


-28.30 


-28.22 


12 C 


-92.16 


-48.68 


14 C 


-105.28 


-74.49 


16Q 


-127.62 


-113.55 


40 Ca 


-342.05 


-340.75 


48 Ca 


-416.00 


-362.82 


208p b 


-1636.45 


-2059.55 



Table 2: Binding energies, in MeV, obtained with the interaction obtained by adding to 
the Bl interaction the tensor terms of ref. M (JSB) and that of ref. |8| (A). The column 
labeled T6 shows the results obtained with the tensor force adjusted to reproduce the 
experimental value of the 12 C binding energy. 





JSB 


A 


T6 


4 He 


-28.22 


-28.22 


-28.22 


12 C 


-47.17 


-47.46 


-92.36 


14 C 


-74.56 


-74.59 


-90.45 


16Q 


-114.10 


-114.10 


-113.53 


40 Ca 


-342.40 


-342.40 


-340.75 


48 Ca 


-364.01 


-363.98 


-377.01 


208 Pb 


-2074.51 


-2074.00 


-2063.93 



Table 3: Single particle energies, in MeV, for 12 C for the interactions used. Since these 
calculations have been done without the Coulomb interaction, proton and neutron single 
particle energies are the same. 



12 C 


Bl 


JSB 


A 


T6 


lSl/2 


-36.81 


-36.66 


-36.55 


-42.02 


1P3/2 


-11.82 


-11.29 


-11.36 


-25.82 
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Table 4: Single particle energies, in MeV, for 14 C. We have indicated with 7r the proton 
and with v the neutron levels. 



14 C 




Bl 


JSB 


A 


T6 


lSl/2 


7T 


-44.83 


-44.93 


-44.93 


-44.24 


1P3/2 


7T 


-18.09 


-17.98 


-18.00 


-25.65 


lSl/2 


V 


-40.55 


-40.62 


-40.61 


-42.21 


1P3/2 


V 


-15.05 


-14.76 


-14.81 


-25.26 


lPl/2 


V 


-15.05 


-15.77 


-15.66 


-16.91 



Table 5: Same as table |3| for 16 0. 



16Q 


Bl 


JSB 


A 


T6 


lSl/2 
1P3/2 
lPl/2 


-49.40 
-21.71 
-21.71 


-49.55 
-21.79 
-21.79 


-49.55 
-21.79 
-21.79 


-49.39 
-21.71 
-21.71 



Table 6: Same as table [I] for 40 Ca. 



40 Ca 


Bl 


JSB 


A 


T6 


lSl/2 


-71.34 


-71.57 


-71.57 


-71.34 


1P3/2 


-45.24 


-45.39 


-45.39 


-45.24 


lPl/2 


-45.24 


-45.39 


-45.39 


-45.24 


ld 5 /2 


-21.19 


-22.04 


-22.04 


-21.95 


ld 3 / 2 


-21.19 


-22.04 


-22.04 


-21.95 


2si /2 


-20.17 


-20.26 


-20.26 


-20.17 
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Table 7: Same as table |]for 48 Ca. 



48 Ca 


Bl 


JSB 


A 


T6 


lSl/2 


TT 


-73.87 


-74.08 


-74.07 


-74.43 


1P3/2 


7T 


-50.81 


-50.56 


-50.53 


-51.76 


lPl/2 


TT 


-50.81 


-51.03 


-51.06 


-51.04 


lds/2 


7T 


-28.95 


-28.72 


-28.70 


-34.07 


ld 3 /2 


TT 


-28.95 


-29.50 


-29.50 


-22.45 


2si /2 


TT 


-26.96 


-27.54 


-27.48 


-26.96 


lSl/2 


;/ 


-69.57 


-69.78 


-69.77 


-70.06 


1P3/2 




-45.44 


-45.49 


-45.47 


-44.96 


lPl/2 




-45.98 


-46.27 


-46.29 


-48.87 


ld 5 / 2 


V 


-23.64 


-23.57 


-23.55 


-25.87 


lda/2 


V 


-24.47 


-24.78 


-24.79 


-21.74 


2si /2 


V 


-22.94 


-23.03 


-23.29 


-22.88 


2f7/2 


V 


-4.19 


-4.06 


-4.05 


-7.57 
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Figure Captions 



Figure 1: Tensor isospin terms added to the Bl interaction. In the panel (a) the curve 
labeled with A represents the term of the Argonne Vu potential ||, that labeled as JSB 
the term of the Jiilich-Stony Brook interaction of ref. ||. The panel (b) shows the 
interaction adjusted to reproduce the 12 C experimental binding energy. 



Figure 2: Proton density distributions obtained with the Bl interaction (full lines) with 
the A interaction (dashed lines), with the JSB interaction (dotted-dashed lines and with 
the T6 interaction (doubly dotted-dashed lines). On the scale of the figure the A and 
JSB results are overlapped. 
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